Optical properties of Fano-resonant metallic metasurfaces on a substrate 



S. Hossein Mousavi, Alexander B. Khanikaev, and Gennady Shvetifl 
Department of Physics and Institute for Fusion Studies, University of Texas at Austin, 
One University Station C1500, Austin, Texas 78712, USA 

Three different periodic optical metasurfaces exhibiting Fano resonances are studied in mid-IR 
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provide a simple description of the observed behavior. We demonstrate that the presence of a 
substrate substantially alters the collective response of the structures suppressing Wood's anomalies 
and spatial dispersion of the resonances. Different types of Fano resonances are found to be affected 
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the validity of the widely used effective medium approaches is re-examined. 
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I. INTRODUCTION 

Optical nanostructures exhibiting Fano resonances 
have attracted a significant attention of the research 
community in recent yearsi~i£. As classical analogues 
of quantum Fano system^, these systems represent el- 
egant tabletop tools for testing fundamental principles 
of physics^. From the applications' point of view, the 
possibility to trap, enhance, and manipulate light in op- 
tical Fano nanostructures and metamaterials is also very 
promising and has already been demonstrated to be ben- 
eficial in sensing and bio-sensing 15 18 , photovoltaics and 
thermo-photovoltaicsi^r— , and slow-light generation 3 . 

Fano resonances in optical metamaterials originate 
from electromagnetic interactions between their con- 
stituents. According to the type of interaction, Fano 
resonances can be classified into two groups: (a) coher- 
ent Fano resonances which are the result of interferences 
between all of the meta-atoms forming a periodic ar- 
ray, and (b) local Fano resonances which originate from 
the complex local structure of individual meta-atoms. 
While for the first class of these resonances, periodic ar- 
rangement of a large number of meta-molecules is cru- 
cial, for the second class, even a single meta-molecule 
may exhibit a Fano resonance 2 -. However, in both cases 
the response of the meta-surface will depend on the di- 
electric environment. For coherent metamaterials, it is 
well-known that the presence of a substrate modifies the 
far-field interaction s 22 ! 23 , which substantially alters the 
collective response. However, even in the case of lo- 
cal Fano resonances, the presence of a substrate can 
significantly modify the interaction between the meta- 
molecule's constituents 24 . Thus, understanding of sub- 
strate effects in both cases is critical for basic under- 
standing of Fano resonances in experimentally relevant 
systems. 

The geometry of the systems with Fano resonances can 
vary from simple designs such as metal nanoparticles 2 
and perforated films 25 27 to structures with complex ge- 
ometries such as oligomers and dolmens^ " 1 1 1 14 . 



The simplest approach to describe the Fano resonances 
is a mechanical toy model of coupled harmonic oscillators 
— &L2£. In many cases, such model provides sufficient in- 
sight into the physics of the system; it allows one to fit 
optical spectra and associate spectral features of the sys- 
tem with the known properties of the mechanical Fano 
model. Nevertheless, such simplified representation of a 
complex physical system neglects various effects relevant 
for real structures. Among these are substrate effects, 
effects of periodicity resulting in a collective character of 
the electromagnetic response of meta-materials, geome- 
try/proximity effects modifying the interactions between 
meta-molecules. 

Another analytical description applicable to some rela- 
tively simple systems, such as arrays of point-like dipoles, 
is offered by a dipole model 5 -* 2 ^"— . Indeed, a few of the 
above-mentioned phenomena, including collective effects 
and Fano resonances, can be adequately described by this 
model 33 . However, there are situations when the dipole 
model is inadequate. For instance, as soon as the scat- 
terers or the metamolecules forming the array become 
comparable in size to the wavelength of light, the dipole- 
dipole interaction mechanism breaks down because in- 
teractions through higher multipoles come into play. In- 
corporating substrate effects also represents a challenge 
for the dipole model. To resolve these problems, the 
dipole model has been extended by including effects of 
higher multipoles 34 , found with the use of numerically 
calculated polarizabilities of metamolecules and by ap- 
plying the scattering-matrix technique 35 ^ 3 - 7 .. However, af- 
ter such generalizations, the dipole model can no longer 
be considered analytical. To fully account for all the 
above-mentioned effects, the approach which is being 
widely used is to apply a powerful yet time consum- 
ing ab-initio numerical solvers of the Maxwell equations, 
or various generalizations of hybrid numerical/analytical 
approaches 3 -^. 

However, there is another very promising semi- 
analytical yet rigorous approach based on the modal- 
matching technique (MMT). MMT has been proven to be 
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FIG. 1. (Color online.) Schematic unit cell of the three struc- 
tures considered in this paper: (a) single-antenna, (b) double- 
antenna, and (c) dolmen metasurfaces. In (c), two unit cells 
(separated by a dashed line) are shown, (d) shows the side 
view of the metasurfaces cladded by the superstrate and the 
substrate. 

very fruitful for describing various electromagnetic sys- 
tems ranging from frequency-selective surfaces and an- 
tenna arrays in radio frequency domain 39 to perforated 
metallic structures 4 ^"— . In both cases, the approach is 
relying on the expansion of a system-specific polarization 
(e.g., the total current in the antennas or the fields inside 
the holes) in the form of a superposition of the eigen- 
modes satisfying appropriate boundary conditions. It is 
important that, because the basis used for the expan- 
sion takes into account met a- molecules' geometry, this 
approach can account for the geometry-specific effects, 
including shape resonances and proximity effects 44,45 . 

In this paper, the semi-analytical MMT based on the 
current expansion— is applied to study Fano resonances 
in two-dimensional periodic arrays of metallic antennas 
on a substrate (Fig. [1]). The emphasis is on arrays of 
complex antennas resonant in the mid-IR part of opti- 
cal spectrum, which is important for bio-sensin g 16 i 18 and 
thermo-photovoltaic applications^"— . MMT is shown 
to capture the meta-molecule geometry as well as the 
substrate effects and is easily expandable to the cases of 
more complex geometries composed of several spatially 
extended scatterers per unit cell. We develop a minimal 



model of the MMT by truncating the electric-current ba- 
sis of the full model and demonstrate that even the result- 
ing analytical model is capable to quantitatively describe 
the optical properties of the structures. The main em- 
phasis of the paper is on the effect of the substrate. It 
is demonstrated that with the model in hand the MMT 
can predict the spectral positions of the resonances and 
describe the interference among the different scattering 
pathways provided by the different metamolecule's reso- 
nances. 

The rest of the paper is organized as follows. In Sec- 
tion II a general MMT based on the Rayleigh and current 
expansions is derived. In Section III a minimal model is 
developed and applied to study a periodic single- antenna 
metasurface (SAM), shown in Fig. [TJi. In Section IV we 
analyze periodic double-antenna meta-surfaces (DAMs) 
[Fig. [T^], which are known to exhibit a Fano resonance 
due to the interference of the modes corresponding to 
symmetric and anti-symmetric charge distributions in the 
antenna pairs. Finally, in Section V a dolmen struc- 
ture formed by three antennas (Fig. []J) and known to 
exhibit a plasmonic analogue of electromagnetically in- 
duced transparency (EIT) 6,7 is studied. 

II. MODAL MATCHING TECHNIQUE: 
THEORETICAL FORMALISM 

In this section the electric current expansion technique 
is generalized to the case of infinitesimally thin metallic 
antennas of finite surface conductivity. Following Ref^, 
we rely on the current expansion, but require the current 
in the plasmonic antennas to be defined by the high- 
frequency conductivity of the metal and the electric field 
right on the antennas' surface. To make the expressions 
more compact and keep the formulation more general, 
the Dirac notation is used. Using the periodicity of the 
structures the tangential components of the fields are ex- 
panded in the superstrate (I) and the substrate (II) in 
the plane- wave basis, 
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where |n,r) = \n Xl n yi r) represents the in-plane (tan- 
gential) electric field of the n-th diffracted plane wave 
with an in-plane wave-number k" = (fc£,/c"), a polar- 
ization state r (s or p polarization) and a wave admit- 
tance in the superstrate (substrate) (Y^ 1 ). The 
wave admittance of an s-polarized plane wave is Y n s = 
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fc"/(Zofco) and that of a p-polarized plane wave is Y n p = 
€ r ko/(Zok™). In the above expressions, k n = 2im x /P x1 

ky = 27rn y /Py, and k n — yj 1 e r k1 — (k n ) 2 are respectively 

the x, I/, and z components of the wavevector k. fco is the 
wavenumber of the wave in vacuum and e r is the relative 
permittivity of the medium in which the wave is propa- 
gating. Z « 376.73 Q is the impedance of vacuum. The 
spatial representation of | n, r) is given by: 

, , exp (ik n x + ik n y) 

( r ni^) = , n /pp v - m*- m, (2) 

V r x±y 

. . exp (i/c^x + ifc n v) 

<r„| n ,p>= + W (3) 

^|| V -Lx-Ly 

The surface current in the plasmonic antennas is ex- 
panded with the use of the basis functions | j™) satisfying 
the boundary condition of vanishing current on the an- 
tenna edges. The electric current in the ra-th antenna 
within the unit cell takes the form 



= £C|ja> 



(4) 



where c™ is the amplitude of the a-th current mode. The 
precise form of the functions | j™) depends on the antenna 
geometry, which will be taken rectangular throughout the 
paper. For this case, the explicit expressions for | j™) are 
given elsewhere.— 

The finite conductivity of the antennas is incorporated 
into the model using the relation 



(5) 



where a = ihe m ko is the effective surface conductivity of 
the antennas of thickness h and permittivity e m 46 . The 
fields in the substrate and the superstrate should satisfy 
the standard continuity boundary condition over the in- 
terface between them 
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Over the antennas, the electric field is continuous and 
Eq. (|6aj) still holds, however a discontinuity in the mag- 
netic field appears due to the finite current in the anten- 
nas 
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Combining these constraints together with the modal ex- 
pansions and using the orthonormality of the basis func- 
tions, we obtain a system of three linear equations: 

C(V - r n ,r) - On,r = £ (n, r IOC (7b) 
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Eqs. ([TK-IZt>) are obtained by multiplying Eqs. flSk-c) 
(from left) by (n, r| and Eq. ([Tfc) is obtained by multi- 
plying Eq. (|5]) (from left) by The scalar product 
(n, r|j) is defined as 
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where j x (x,y) and j y (x,y) are the x and components 
of the current profile |j). Eqs. jTMZh) can be solved to 
obtain the amplitude of the scattered fields tn, T and r n ^ r 
in terms of the antenna current modes c™ 
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where the amplitude of the a-th current mode in the m- 
th antenna c™ is 
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In Eq. ([T2j) ? S™™/ is the Green's function that describes 
the cross-talk between the a-th current mode of the ra- 
th antenna with the at-th current mode of the m/-th an- 
tenna in the unit cell. This coupling is mediated by the 
plane waves of all (propagating and evanescent) diffrac- 
tion orders n as given by Eq. (|13aj) . %™ is the direct 
coupling strength of the a-th current mode of the m-th 
antenna to the external field. 
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In addition to the scattering characteristics, the eigen- 
modes of the structure can also be determined by solving 
the secular equation 



det{a 1 S aia ,S rt 



at i — U « 



(14) 



In general, Eq. ([H]) can be satisfied only at complex 
frequencies uo — uo r + iuui, showing that the eigenmodes 
have a finite lifetime due to either Ohmic losses (when 
5fte{cr} ^ 0) or radiative deca y 42 1 43 . A real eigenvalue 
uj may be ideally achieved only in the limit of lossless 
meta-surfaces and for fen > fc, which ensures absence of 
absorption and radiation. While in this paper we assume 
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a two-dimensional periodic array of the antennas, but the 
results can be generalized to the case of an individual an- 
tenna. In that case, the sum over the discrete set of plane 
waves in the above expressions, should be replaced 
by an integral over the continuum J J dk x dk y . A7 



III. COLLECTIVE RESPONSE OF 
SINGLE- ANTENNA METASURFACES 

While Eqs. ()10til2|) fully characterize the optical prop- 
erties of the system, they have to be solved numerically. 
However, in the case when the wavelength of the inci- 
dent light matches the resonant wavelength of a partic- 
ular current eigenmode truncation of the current 
modes' basis can provide simple and instructive expres- 
sions 48 . Indeed, it has been shown for the case of per- 
forated metal films that near the cut-off frequency of a 
particular waveguide mode such a "minimal" model can 
provide a very good approximatio n 41 ^ 9 . As will be shown 
later this situation holds for the antenna geometries in 
the frequency range studied here. It is worthwhile men- 
tioning here, however, that as dimensions of antennas 
decrease and they become increasingly subwavelength, 
the convergence of the MMT deteriorates and more and 
more current modes should be considered. 

By limiting the current basis to one mode per antenna 

I J i 1 ) = ^/? cos (cP^)^ we assume a dominant role for 
the fundamental (dipolar) antenna-mode. Such a trun- 
cation allows us to get simple analytical expressions for 
the transmission and reflection amplitudes which quali- 
tatively explain the system's response yet quantitatively 
match the numerical results. 

First we consider a structure with a single antenna per 
unit cell with large length-to-width aspect ratio so that 
the fundamental current mode | j) = | jj;) with the current 
along the long antenna dimension dominates. In this 
case Eq. ([T2]) can be analytically solved for the current 
amplitude 
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where X i(k|| , u) = 2Y^ r /(Y^ r + F n r ) (j 1 0, r ) is the cou- 
pling efficiency of the antenna to the external field, the 
sum 5ii(k||,u;) = SH(k\\,uj) as defined in Eq. (j!3aj) . 
^ext = ^o,r and the rest of the z n ,r's are assumed to be 
zero. In section III and IV we assume that the incident 
polarization is along the long dimension of the antennas. 
In the minimal model, Su is explicitly given by: 
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FIG. 2. (Color online.) (a) Zeroth-order normal-incidence 
transmission (solid lines) and amplitudes of the fundamental 
current mode (dashed lines) for the single- antenna metasur- 
face for the symmetric (ei = en, black lines) and asymmetric 
(ei 7^ en, red lines) claddings, (b) Zeroth-order transmission 
of the single- antenna metasurface at the Wood's anomaly as 
a function of the dielectric contrast ei — en (solid line) and 
the bare (without antennas) interface transmission (dashed 
line). The structure parameters are as follows: a x = 0.3 /xra, 
a y — 1.5 /im, P x — 1.8 /im, P y — 1.8 firn. and en — 12. 



where the first term in the square bracket corresponds 
to the s-polarized diffraction orders and the second 
term corresponds to the p-polarized diffraction orders. 
sinc(t) ee sin(t)/£ can be approximated by unity for thin 
antennas (t = k x a x /2 w 0). We emphasize that this sum 
converges much faster compared to that in the dipole 
model. In the framework of the minimal model the trans- 
mission and reflection coefficients of the single- antenna 
meta-surface assume a simple form: 
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Eqs. (fT7( [T8 ]) explicitly show two scattering pathways: (i) 
the direct transmission (reflection) through the dielectric 
interface between the substrate and the superstrate with- 
out interacting with the antennas and (ii) the field radi- 
ated by the antennas due to the excitation of the current 
mode. These two scattering pathways and the presence 
of the Wood's anomal y 50 ! 51 lead to a Fano interference 
resulting in an asymmetric lineshape typically observed 
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in antenna arrays (Fig. [2^) and systems with extraordi- 
nary optical transmission. The condition of the vanishing 
denominator in the Eq. (fT5]h ct _1 (cj) + 5n(k||,a;) = 0, 
approximately describes the eigenmodes dispersion of the 
antenna array (cf. Eq. ([T4]) for the exact expression). 

0.1 \ 1 
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FIG. 3. (Color online.) Direct transmission at normal in- 
cidence calculated using three methods: (i) first-principles 
COMSOL simulations (circles), (ii) full MMT model (solid 
line), and (iii) the minimal MMT model (dashed line). Struc- 
ture parameters: same as in Fig. [2] 

Note that Eq. ([T5]) strongly resembles that of the ef- 
fective polarizability of an array of dipoles in the dipole 
model^. In our case, ^m{l/a(uj)} plays the role of 
"plasmonic polarizability" of the antennas while the role 
of 5ii(k||,co>) is similar to that of the dipole lattice sum. 
Note that because the modal-matching technique is based 
on current expansion, the polarizability and the lattice 
sum have tt/2 phase shift as compared to the dipole 
model. It can be shown that when the dipole sum of 
the dipole model is written in the reciprocal space, the 
sums of both models have similar terms diverging at the 
Wood's anomalies. However, in contrast to the dipole 
model, the antenna model fully considers the substrate, 
antenna shape and finite conductivity. Here we limit our 
consideration to the mid-IR domain where antennas can 
be considered as perfectly conducting (1/<t=0). Effect 
of the finite conductivity of the antennas on Fano reso- 
nances will be considered elsewhere. 

Figure [3] compares the results obtained with the use 
of the minimal model, the full MMT with a sufficient 
number of current modes (such that the convergence is 
reached), and the full- wave COMSOL Multiphysics sim- 
ulations. The minimal model shows very good agree- 
ment with both of the exact techniques. The inclusion 
of the higher-order modes only insignificantly changes 
the spectra. Some discrepancy between the full MMT 
and COMSOL results is observed in the long wavelength 
range which is attributed to the rounding of the antenna 
corners that was done in COMSOL to minimize the nu- 
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FIG. 4. (Color online.) Angular-resolved zeroth-order (di- 
rect) transmission through the single-antenna metasurface in 
(a) symmetric (ei = en) and (b) asymmetric (ei = 1, en = 12) 
cladding. The parameters of the structure are the same as in 
Figs. [2]and El 

merical error due to spurious singularities at sharp cor- 
ners. 



A. Substrate effects on Wood's anomalies 

The substrate effect appears due to the finite dielectric 
contrast Ae = e\ — en between the superstrate and sub- 
strate claddings, and in the MMT it is described by the 
wave admittances (Y njT ) entering the expression for the 
sum Sn [Eq. (|13aj) ]. Predictions of the minimal model 
and the exact calculations clearly show that the presence 
of a substrate strongly affects the scattering character- 
istics of the antenna array (Fig. [2]). This effect can be 
fully understood from the analytical expressions of the 
minimal model. Let's first consider the case of a Wood's 
anomaly that indicates the onset of a diffraction order 
and appears in the transmission spectrum as a maxi- 
mum (for symmetric cladding, e\ = en) or as a kink (in 
the general case) and represents the sharpest feature of 



the spectrum. Finite optical contrast between the sub- 
strate and the superstrate makes the sum Sn converge 
at the onset of the s-polarized diffraction orders which 
manifests as the suppression of the Wood's anomalies. 

Indeed, as can be seen from the expression for 
Sii(fc||,a;), in the case of a symmetric cladding (ei=en) 
there is a divergent term corresponding to the onset of 
a s-polarized diffraction order |n, s); when the diffrac- 
tion order experiences a transition between the evanes- 
cent and propagating regimes Y n s ~ k™ vanishes and the 
lattice sum Su(k\\,uj) diverges as l/kf. This gives rise to 
the Wood's anomaly and vanishing of the current mode in 
the antennas c(k\\,(j) —> 0, as illustrated by black dashed 
curve in Fig. [2k- Thus, at the Wood's anomaly, antennas 
are inactive and invisible (to, r = 1 and ro, r = 0) as can 
be seen from Eqs. (|17til8|) . 

In the case of an asymmetric cladding, e\ ^ en, the 
same term of the lattice sum that was divergent for the 
symmetric case, assumes the form l/(k^ 1 + k™ ,u ) and, 
since ^ &"' n , the lattice sum fin (fey, a;) never di- 
verges. As shown in Fig. [2^ by red dashed curve, even at 
the Wood's anomaly, the current in the wires c(ky , ujwa), 
does not vanish and the antenna array scatters the radi- 
ation. 

The fact that antennas remain polarized and con- 
tribute to scattering at the Wood's anomaly for any finite 
value of the optical contrast Ae is illustrated by Fig. [2d, 
where the transmission at the wavelength corresponding 
to Wood's anomaly of the antenna array on the substrate 
(solid line) is plotted along with that of the bare (i.e., no 
antennas) interface (dashed line). The two curves in- 
tersect only for the case of a symmetric cladding and 
the discrepancy can be significant for a large contrast. 
Thus, the peak transmission of the antenna array on the 
substrate never reaches that of the bare interface. It is 
also remarkable that the transmission of the antenna ar- 
ray appears very asymmetric with respect to the sign of 
the contrast and the difference can reach tens of percent. 
This asymmetry is caused by the presence of the Wood's 
anomaly of the second medium whose diffraction order 
dominates in the regime of e\ > en. 




FIG. 5. (Color online.) (a) Comparison between the exact 
treatment of the substrate and the results of the homogenized 
antenna's environment according to the effective medium the- 
ories: (i) e e ff = (ei + en)/2 (dashed-dotted line) and (ii) 
n eff — (ni + nn)/ 2 (dashed line), (b) Deviation of the res- 
onance spectral position (AA/A r ) predicted by the effective 
medium approaches from the exact result A r shown in the 
inset. 



B. Substrate's influence on collective antenna 
resonances 

The spectral position of the collective dipolar antenna 
resonance corresponds to the minimum of the transmis- 
sion or the maximum of the reflection spectra. It can 
also be found from the zeros of the imaginary part of 
the denominator in Eq. [15l In this paper, the length of 
the antennas are such that the observed collective mode 
stems from the A/2 resonance of an individual antenna 
(as can be seen from the charge distribution plotted in 
the inset to Fig. [3j). As can be seen from dashed curves 
in Fig. [2^l, it also corresponds to the maximum of the 
current amplitude c excited in the antennas. Note that a 
good agreement between the MMT and COMSOL spec- 



tra near the resonance in Fig. [3] indicates that our model 
fully accounts for the modified polarizability of the an- 
tennas caused by the presence of the substrate. 

The presence of the substrate and suppression of the 
Wood's anomalies dramatically affect the angular disper- 
sion of the resonance and change its nature. Figure HJa) 
shows that for the case of a symmetric cladding we ob- 
serve the effect of "dragging" of the mode by the Wood's 
anomal y 32 ! 33 , which results in its strong spatial disper- 
sion. This strongly dispersive character of the mode is 
a manifestation of its collective origin. Indeed, the long- 
range interaction among the individual metamolecules re- 
sults in the appearance of the collective modes which are 
spectrally displaced and sharper as compared to those 
observed in the individual metamolecules 32,33 . However, 
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for the case of an asymmetric cladding, shown in the 
Fig. 0}}, the dragging effect is diminished and the reso- 
nance crosses the Wood's anomalies, every time reducing 
its quality factor due to the opening of additional ra- 
diative channels. The resulting non-dispersive behavior 
and lower quality factor of the mode imply that the reso- 
nance has lost its collective character and represent local 
excitation modified by the complex electromagnetic envi- 
ronment. Note that for the confined mode of the periodic 
metasurfaces, i.e. non-leaky spoof plasmons, presence of 
the substrate and associated suppression of the collective 
behavior may result in the disappearance of the mode 36 . 

By following the resonance one can also track its spec- 
tral position as a function of the dielectric contrast Ae. 
Inset to Fig. [SJd shows the corresponding spectral shift 
of the resonance. The frequency of the resonance shows 
nearly linear dependence on the contrast. We emphasize 
that the effect of the substrate on the antenna polariz- 
ability cannot be rigorously described by a basic dipole 
model, but can be phenomenologically incorporated into 
it. One of the most common approaches to include the 
optical contrast is to use an effective medium approach 52 
and approximate the permittivity or the refractive in- 
dex of the background medium by e e ff = (ei + en)/2 or 
n eff — (ni + nn)/2. The validity of such approxima- 
tions can be tested by comparing the spectral position of 
the resonance calculated for the antennas embedded into 
the effective medium with the exact MMT result. Fig- 
ure [3i shows the deviations of the resonance frequency as 
predicted by the effective medium theories from the ex- 
act result. One can see that the agreement between the 
effective medium approaches and the exact calculations 
is rather good. The homogenization of the permittivity 
e e ff = (ei + €n)/2 is closer to the exact result and pro- 
vides a good estimate even for the case of a large contrast. 
However, it can be seen from Fig. [5^ that the effective 
medium approach fails in predicting the transmission and 
reflection, especially close to the Wood's anomaly. 
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FIG. 6. (Color online.) Comparison of the zeroth-order 
transmission spectra of a double- antenna metasurface on the 
substrate for 50° angle of s-polarized incidence calculated 
with the use of COMSOL Multiphysics, and the minimal and 
full MMT models. The structure parameters are as follows: 
a x — 0.3 /im, a y — 1.5 fim, P x — P y — 1.8 //m, D x — 0.6 /im, 
ei — 1, and en — 12. 



Using the unitary transformation 

Csub = ci exp(ik x D x /2) - c 2 exp(-ik x D x /2), 

c sup = ci exp(ik x D x /2) + c 2 exp(-ik x D x /2), (20) 

the basis of the two fundamental electric currents can be 
changed to the more instructive basis of the sub-radiant 
and super-radiant current modes. In the new basis, the 
system is described by the matrix equation 



Su + A in 
—in Su — A 



°sup 
^sub 



(21) 



IV. FANO RESONANCE IN 
DOUBLE- ANTENNA META- SURFACES 



Next, the MMT for the double- antenna meta-surface 
(DAM) shown in Fig. [Ud is considered. In this case, the 
truncation of the basis to the fundamental modes in both 
antennas, results in a minimal model with a 2x2 matrix 
equation for the amplitudes c\ and 



S\\ S12 

S21 Su 



(19) 



where S mrn , = S 7 ^' represents the effective Green's 
function given by Eq. (fT3a|) . and x=2F I r /(Y I /r + F n T ) 

exp(ik x D x /2) |0,r^ is the coupling strength of the 

incident light (|0,r)) to the fundamental current modes. 



where A=— 1/2 [S12 exp(ik x D x ) + S21 exp(— ik x D x )] and 
n=i /2[S 12 exp(ik x D x ) — S 21 exp(—ik x D x )]. These modes 
are the result of the radiative coupling of the dipolar 
resonances of the antenna pairs^^^. The super-radiant 
and sub-radiant modes have distinct symmetric and anti- 
symmetric charge distributions, while the currents corre- 
sponding to them are collinear and anti-collinear, respec- 
tively, as illustrated by the inset to Fig. [6] The modes 
interaction makes them blue-shifted (super-radiant) and 
red-shifted (sub-radiant) with respect to the electric- 
dipolar resonance of SAM. Another consequence of such 
hybridization and different symmetry of the modes is 
their different radiative coupling efficiency. The sub- 
radiant mode is not directly coupled to the incident light. 
The latter fact is reflected in the mode's name and di- 
rectly follows from the zero radiative coupling strength 
on the RHS of Eq. (j2T|). 

With the use of Eqs. (|10|lip we obtain a set of minimal- 
model expressions for the transmission and reflection co- 
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FIG. 7. (Color online.) (a) Zeroth-order s-polarized transmis- 
sion spectra of the double-antenna meta-surface for the sym- 
metric (ei=en=12, blue line) and asymmetric (ei=l, en=12, 
red line) claddings, (b) Amplitude of the sub-radiant (dashed 
lines) and super-radiant (solid lines) current modes. The pa- 
rameters of the structure and the incidence angle are the same 
as in Fig. [6] 



efficients: 



2Yl 



yl + yll 

n.r 1 n.r 



yl _yll 
n.r n.r . 

yl 7 yll Zr 

n.r ' n.r 



n,r|j (1) )exp(-ifc XJ D x /2) 

yl T yll Csu P ' 

n.r "r n.r 

(22) 

n,T|j«)exp (-2^1^/2) 

C SU p. 



yl + yll 
n.r 1 n.r 



(23) 



Note that because of their coupling, c sup and c su b 
are not the true eigenmodes of the system. The 
true eigenmodes which will be referred to as dipo- 
lar d and quadrupolar q can be obtained using the 
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FIG. 8. (Color online.) Deviation of the spectral position 
(AX/Xr) of the dipolar (red lines with circle markers) and 
quadrupolar (blue lines with triangle markers) resonances pre- 
dicted by the effective medium approaches from the exact re- 
sults X r shown in the inset. 



transformation: 33 



, A + X 
a — — — — — c. 



IK 

2X ~ sup + 2X 
in A + X 



q- 2X c sup 



2X 



Csub? 



(24) 



where X = VA 2 + k 2 . The modal coefficients d and q 
are related to the incident field through 



, 1 + A/X 

d = XEext 



s d 

in/X 



XEe. 



Sg = S u - X, (25) 



In the particular case of normal incidence, the 
quadrupolar mode q exactly coincides with the sub- 
radiant mode c su b since its coupling to the super-radiant 
mode c sup vanishes (k —> as ky —> 0). However, 
at finite angles of incidence, this mode acquires a fi- 
nite electric-dipolar moment and its radiative coupling 
and bandwidth gradually increase. The dipolar mode, 
in contrast, is always strongly radiatively coupled and 
is spectrally broad at any incidence angle. Eqs. (|224I23|) 
describe a Fano resonance of a single "continuum" inter- 
acting with two "discrete" dipolar (d) and quadrupolar 
(q) resonances. However, the bandwidth of the last two 
is very different due to their different radiative coupling. 

The transmission spectrum of DAM for an oblique 
incidence of 50° calculated using the minimal model, 
full MMT and COMSOL solver are plotted alongside in 
Fig.[6l As before, the different techniques show very good 
agreement, which implies the applicability of the analyt- 
ical minimal model even for the description of light scat- 
tering by more complex metasurfaces such as DAM. Two 
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resonances, a broad one due to the excitation of the dipo- 
lar mode and a narrow one due to the excitation of the 
quadrupolar mode, are clearly seen in the spectra. As ex- 
pected, the narrow resonance has a strongly asymmetric 
shape typical for Fano resonances. 

A. Dipolar mode 

Because the dipolar and quadrupolar resonances are 
spectrally separated and the quadrupolar resonance is 
spectrally narrow, first we can consider the dipolar reso- 
nance individually. It has been recently shown that the 
collective dipolar resonance in DAM appears to have spa- 
tial dispersion [o;(k|i)] very different from that found in 
SAM 33 . Using the basic dipole model, this difference 
was shown to be due to the suppression of the long- 
range interactions in DAMs. MMT provides analogous 
results. In the present case, however, the role of the 
substrate and the antenna geometry are fully taken into 
account. In DAM, the Wood's anomalies would have 
corresponded to the divergences (or local maxima for an 
asymmetric cladding) of the effective Green's function 
Sd = Sn + X playing the role of the lattice sum in 
the dipole model or the Green's function Sn of SAM. 
However, it can be shown that exactly at the frequency 
where the Wood's anomaly is expected, the diverging 
terms in Sn and X exactly cancel out eliminating the 
divergence, thereby suppressing the spectral features as- 
sociated with the Wood's anomaly 33 . Thus, in the DAM 
case, Wood's anomaly appears to be suppressed even for 
a symmetric cladding (Ae = 0). Therefore for the asym- 
metric cladding, the Wood's anomaly is suppressed by 
two mechanisms: one due to substr ate/super str ate con- 
trast just as in the case of SAM, and another one due 
to disconnected topology of DAM. Note, however, that 
despite such suppression, one might still observe some 
variations of reflectivity in the vicinity of the expected 
Wood's anomalies, which is a result of the opening of the 
diffraction channel and change of the mode's radiative 
lifetime. 33 . 



B. Quadrupolar mode 

Now we focus on the quadrupolar resonance. The nar- 
row spectral region where the transmission undergoes a 
rapid variation is especially interesting because of the 
strong field enhancement and slow light regime reported 
earlier 3,4 . The asymmetric shape of the spectrum is the 
result of the Fano interference between the quadrupolar 
resonance and a background provided by two transmis- 
sion channels: (i) the direct transmission without inter- 
action with the array, and (ii) the radiation scattered by 
the dipolar resonance. This background is nearly con- 
stant and is featureless over the spectral width of the 
quadrupolar resonance where the transmission changes 
from its maximal value to zero. It can be analytically 



shown that the transmission peak corresponds to the di- 
vergence of c su b and zero of c sup (Fig. 03). Because only 
the super-radiant mode (c sup ) is coupled to the radia- 
tion, the antennas do not radiate at this frequency and 
the transmission acquires a universal value of the trans- 
mission of the bare interface between the substrate and 
the super str ate as seen from Eqs. ()22I23|) . Note that this 
behavior is different from that found for another type of 
Fano resonances, Wood's anomalies in SAM, when for an 
asymmetric cladding configuration, the radiative current 
mode always had some finite amplitude. 

In the previous section, it was demonstrated that the 
effective medium approach accurately predicted the spec- 
tral position of the dipolar resonance for SAM. Here 
we test this approach for both cases of the dipolar and 
quadrupolar resonances in DAM. At first, the inset to 
Fig. [8] shows the spectral position of the resonances cal- 
culated by the full MMT approach as a function of 
the dielectric contrast, and Fig. [8] shows the deviation 
from the exact result. The effective theories match well 
with the exact results for both quadrupolar and dipo- 
lar resonances, while the permittivity homogenization 
e e ff = (ei + en)/2 is more accurate than the index ho- 
mogenization n e ff = (ni + Tin)/ 2. However, just as in 
the case of SAM, none of the homogenization approaches 
succeed in predicting the transmission/reflection spectra. 

Finally, we observe that in contrast to the dipolar 
mode, the quadrupolar mode in DAM is strongly affected 
by the Wood's anomaly. From the expression for its am- 
plitude Eq. ([25]) one can see that for the quadrupolar 
mode the divergence in the Green's function Sn — X does 
take place in the case of a symmetric cladding ei = en. 
The effect of the Wood's anomaly on the quadrupolar 
resonance of the DAM is expected to be especially signif- 
icant for large angles of incidence, when they approach 
each other. This results in a strong spatial dispersion 
of the quadrupolar mode observed earlier both in mid- 
and near-IR plasmonic DAMs 33 . For the case of an 
asymmetric cladding the effect of the Wood's anomaly 
on the quadrupolar resonance is again suppressed due to 
the mismatch of the wave admittances in the substrate 
and the superstrate and all the arguments about the sup- 
pression of the Wood's anomalies used in the previous 
section for SAM are applicable here. 



V. ELECTROMAGNETIC ALLY INDUCED 
TRANSPARENCY 

In this section a periodic dolmen structure with three 
antennas in its unit cell 6 ^^*^ [shown in Fig. QJ] is 
studied. This design was first introduced to mimic EIT 
in plasmonic structures for normally incident x-polarized 
light 6 . Two vertical antennas on the dolmen structure 
are responsible for formation of the quadrupolar mode 
analogous to that described in the previous section. This 
mode plays the role of the sub-radiant mode of EIT. The 
horizontal antenna provides a spectrally broad dipolar 
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response. Its length is chosen in such a way that the 
frequency of the dipolar resonance is matched to the 
frequency of the sub-radiant (quadrupolar) mode. The 
position of the horizontal antenna within the unit cell 
defines the degree of symmetry breaking of the meta- 
molecule and the intensity of coupling between the two 
modes. The sub-radiant mode is completely decou- 
pled from x-polarized incident radiation if the antenna 
is placed equidistantly between the two nearest vertical 
double antennas, D y = 0. However, as soon as the hor- 
izontal antenna is displaced vertically from this position 
(D y ^ 0), the sub-radiant mode couples to the incident 
radiation indirectly through the dipolar mode of the hor- 
izontal antenna. This coupling to the sub-radiant mode 
results in a transmission/reflection spectrum typical for 
the EIT: a sharp transmission resonance embedded into 
a region of highly reflecting background provided by the 
dipolar antenna mode. Note that for the polarization 
under study the dipolar mode of two- vertical antennas is 
not excited since there is no ^-component of the incident 
electric field. In principle, this mode could be excited in- 
directly provided that the symmetry of the metamolecule 
is reduced further by displacing the horizontal antenna 
in the x-direction 3 . However, this situation will not be 
considered here. 

Within the minimal model, limiting the current ba- 
sis to only fundamental current mode for each antenna 
and neglecting finite conductivity of the metallic anten- 
nas (1/cr = 0), the EIT structure can be described by a 
3x3 matrix equation relating the current amplitudes to 
the incident electric field: 



(26) 



where S mrn / = S^ m/ represents the effective Green's 
function given by Eq. (|13ap . and x = < J3 10 > 
2F I r /(F I r + Yq 1 ) is the coupling strength between the 
incident light and the horizontal antenna (labeled as "3"). 
The RHS of Eq. ([26|) reflects the fact that the vertical an- 
tennas are not coupled to the x-polarized incident light. 
On the other hand, the presence of off-diagonal matrix 
elements in the LHS of this equation indicates that they 
are coupled to the horizontal antenna. It follows from 
the symmetry of the structure that the basis of three 
currents is redundant for the case of x-polarized normal 
incidence, under which the dipolar mode of the antenna 
pair D y cannot be excited. By performing the unitary 
transformation D y = c\ + C2, Q = c\ — C2 and D x = C3, 
the non-interacting D y mode can be eliminated. The 
truncated governing equation for the super-radiant mode 
D x and the sub-radiant mode Q assumes the form: 



Su 


Sl2 


Sl3 " 




Cl 




"O" 




Su 


-5ia 




C2 


= X^ext 





S\3 


— Sl3 


£33 




_C3_ 




1 



Sq k 

K S Dx 



D x 



XEe 



(27) 



where Sq = l/2(5n - S12), Sd x = £33, and the coupling 
between the modes is given by k = Sis. Note, that unless 




FIG. 9. (Color online.) Substrate effect on the normal- 
incidence zeroth-order transmission spectra of the dolmen 
metasurface for different values of symmetry-breaking param- 
eter D y . Horizontal antenna's dimensions are a x = 1.79 /xm, 
a y = 0.6 /im. Vertical antennas are identical with a x — 0.6 /im 
and a y = 1.7/xm and separated by D x = 1.2 /xm. ei = 1, 
en — 12, and P x —P y —A /xm. 



the horizontal antenna is placed symmetrically (in which 
case k = 0), the modes D x and Q are not eigenmodes of 
the EIT structure, since they are coupled to each other 
through the off-diagonal matrix elements in the LHS of 
Eq. ([27]) . This coupling gives rise to hybridization of the 
sub-radiant mode Q and the super-radiant mode D x and 
formation of new mixed states - quasi-dipolar D x and 
quasi-quadrupolar Q. Both of these modes are coupled 
to the incident electromagnetic field, but have very dis- 
parate coupling efficiency. The true eigenmodes of the 
structure can be found through an eigen-decomposition 
procedure, D x = 1/(1 + r] 2 ) D x + 77/(1 + rj 2 ) Q and Q = 
l/(l+77 2 ) Q-v/(lW) D x where rj = (S Q -S Dx - A)/2« 
[that tends to zero in the case of a vanishing coupling 
k -> 0] and A = y^(S Q - S Dx ) 2 +4^ 2 . Then the equa- 
tions for the eigenmodes' amplitudes assume the form 



D x 



2/(l + r7 2 ) 
S Q + S Dx -A 

-2r7/(l + 77 2 ) 



X^ext, 
X^ext, 



(28) 
(29) 



and the reflection/transmission coefficients of the struc- 
ture can be found from the following equations: 



Yl T -Y£ T , (n,r|j 3 ) 



yl 1 yll y 1 +Y 11 
± n,r 1 ± n,r ± n,r 1 ± n,r 

, <n, | j 3 > 



2Y1 



y 1 +r n y 1 +r n 

n.r 1 n,r n,r 1 n.r 

D x = D x - rjQ. 



D x , 



D x , 



(30) 

(31) 
(32) 
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The transmission spectra of the structure are shown in 
Fig. [9j The different curves demonstrate how the radia- 
tive coupling of the sub-radiant mode increases when the 
symmetry of the metamolecule is gradually reduced by 
increasing the parameter D y . This manifests, at first, as 
an appearance of the sharp and narrow EIT transmission 
peak and a gradual increase of the EIT peak bandwidth. 
It is remarkable that regardless of the symmetry breaking 
degree, the maximum of the EIT peak always tends to 
the same value defined by the transmission of the bare di- 
electric interface between the substrate and superstrate. 
Therefore we can claim that this is rather general rule for 
all Fano resonant systems, where a Fano resonance orig- 
inates from the structure of the metamolecule and not 
from Wood's anomaly as in the case of SAM. 

VI. CONCLUSIONS 

We investigated the effects of a substrate on three dif- 
ferent mid-IR Fano metasurfaces and found that a fi- 
nite refraction-index contrast between the top and bot- 
tom claddings dramatically changes the optical response 
of the structures. In addition to the expected spec- 
tral shift of the resonances induced by the substrate, a 
dramatic change in the collective behavior of the sys- 
tems was found. The widely used effective medium 
techniques were tested. While satisfactorily predict- 
ing the spectral positions of the resonances, they failed 
to describe other scattering characteristics of the struc- 
tures. The Wood's anomalies and modes' dispersion were 
found to be strongly affected by the presence of a sub- 
strate. Divergences corresponding to the onset of the 
s-polarized diffraction orders disappear whenever the re- 
fractive indices of the substrate and the superstrate are 
mismatched, resulting in the suppression of the spectral 
features associated with the Wood's anomalies and giv- 



ing rise to anomalously fiat dispersion of the modes. As 
the Wood's anomaly and the surface resonances of the 
periodic metasurfaces can be considered as a Fano reso- 
nance due to the collective interaction of metamolecules, 
it can be concluded that the presence of a substrate de- 
stroys the Fano picture. As a consequence, the dramatic 
effect of the onset of the diffraction orders reduces to the 
less dramatic change of the modes radiative lifetime. The 
nature of the modes also changes as they lose their collec- 
tive character and start resembling local resonances yet 
modified by their complex environment. 

In addition to this collective Fano resonance, we also 
studied the effect of the substrate on Fano resonances 
originating from the local geometry of the metamolecules. 
Two types of Fano systems were considered with a sharp 
quadrupolar resonance detuned from (double- antenna 
metasurface) and tuned to (dolmen metasurface) the 
dipolar resonance providing the Fano background. In 
both systems the Fano resonances survived after the in- 
troduction of the substrate but were affected due to the 
modification of the frequencies of both low-Q (dipolar) 
and high-Q (quadrupolar) resonances as well as addi- 
tional reflection from the substrate. A universal behavior 
was discovered for both systems: (a) the peak transmis- 
sion reached the value corresponding to the bare inter- 
face, and (b) antennas became effectively invisible at the 
frequency of the Fano resonance. These findings will be 
useful for designing non-dispersive photonic devices with 
a wide-angle optical response. 
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